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Abstract 

Texture advection is an effective tool for animating and investigating 2D flows. In this paper, we 
discuss how this technique can be extended to 3D flows. In particular, we examine the use of 3D and 4D 
textures on 3D synthetic and computational fluid dynamics flow fields. 
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1 INTRODUCTION 

This paper provides some initial results of using texture advection techniques for visualizing 3D flow 
fields. Animation is an indispensable tool for understanding the dynamical properties of velocity fields. 
One has several options on what and how to animate the field. For example, one can animate arrow 
plots and they have some limited effectiveness for 2D fields. But they quickly become too cluttered for 
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3D velocity fields. Alternatively, one can trace and animate particles in 2D and 3D, and for both steady 
and unsteady flow fields [4, 9], While advecting particles are effective means of animating and showing 
the dynamics of the field, the particles are discrete and pointilistic. There tire techniques that attempt to 
provide a surface-like property to these discrete particles [10, 12], Rather than starting with particles, 
our approach in this paper is to investigate other representations. 

Our goal is to produce an animation of the 3D flow field in such a manner that one might observe with 
cloud patterns e.g. large mass of swirling clouds in cyclonic events, or how the jet trail dissipates as it is 
affected by the wind, etc. As such, we look at ways of representing 3D cloud like textures that we can 
then let the flow field advect to generate our animation. In our initial investigation, we use simple fuzzy 
spheres as our “clouds” to see if the idea merits investigation. We also report on how these spherical 
cloud textures advect and deform over time using both a synthetic and a computational fluid dynamics 
(CFD) data set. 

In the next section, we provide an overview of texture advection as a mechanism for 2D flow visual- 
ization. In particular, how one can work around some difficulties such as providing continuous influx of 
new textures at inflow boundaries, or source critical points. We then discuss how we generated our 3D 
and 4D spherical cloud textures. In section 4, we discuss how the advected textures are volume rendered. 
Results of this technique on the two data sets are presented in Section 5. We finally summarize the work 
and discuss our future directions. 

2 2D TEXTURE ADVECTION 

Given a texture image and some underlying 2D flow field, the idea is to successively distort the input 
image over time such that it reveals the flow behavior when animated. We consider each pixel in the 
image as a massless particle and compute its path in the flow field using a second-order Runge-Kutta 
integration method. Suppose that there are n time steps in the animation. Then, at time step fj, where 
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i < n, each pixel is integrated backward by ti time steps. We then replace the current texture at each 
pixel by the texture value of the pixel where the integration ends and save the distorted texture image 
at time t t . To animate the flow, we simply playback the saved images from all n time steps. 

This texture advection method is a simple and effective approach to animate the 2D flow. However, 
the method may not work well for some flows. Specifically, there are two cases where the method may 
fail to depict the flow behavior accurately: (a) if the flow has source critical points, where the flow 
appears to be repelling away from the critical point, and (b) if there are incoming flows from the grid 
boundary. In our earlier work [5], we have proposed several methods to deal with these two cases. Here 
we briefly describe the problems and give an overview of our approaches. The problem described here 
can occur in 3D flows as well. 

In the case where the flow has one or more source critical points, some particles may reach the one of 
the critical points, where the flow velocity is zero. The integration algorithm would then end prematurely 
because the velocity is zero. The question then is what texture value should be used? If we simply use 
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whatever texture is at that critical point, then that texture value would expand over time in the animation. 
In the case of a source critical point, it is very likely that particles which are at an equal distance away 
from the critical point will reach the critical point at the same time. The textures on these particles 
would then form a circular pattern. Hence, we would see the a constant texture expanding radially from 
the critical point. This is undesirable because the animation then would not reveal the flow direction 
accurately. 

In the case where there are incoming flows from the boundary, some particles will reach the boundary 
during the integration. Analogously, what is the texture value to use? If we use whatever the texture is at 
the boundary, then we would see constant bands of texture coming from the boundary. This is also not a 
desirable outcome. 

We have proposed two methods to deal with these two situations. The first method is based on 3D 
procedural textures During the integration, when the particle reaches either the boundary or the source 
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critical point, where the integration ends prematurely, we step into a 3D texture volume by some depth 
to determine the new texture value. The depth is based on the number of time steps remaining when 
the particle integration ends and the texture color is a function of this depth value. We define the 3D 
procedural texture such that the color of the textures changes as one steps into the 3D texture volume. 
This provides an effect of dynamic continuously feeding texture from the boundary as well as source 
critical points. 

Another method that we have proposed is streamline cycling [5] For a subset of pixels in the image, 
we consider each pixel as the seed of a streamline, where we integrate the particle forward and backward 
until it (a) reaches the boundary, (b) reaches a critical point, or (c) has exceeded the maximum length of a 
streamline. The latter case is necessary to prevent very long streamlines in circular flow regions. During 
the streamline computation, we save all the textures at each integration step. To animate the flow, we 
simply shift the textures along the streamline by a fixed amount at each time step of the animation. We 
found that it is not necessary to compute one streamline per pixel. It is sufficient to seed the streamlines 
every two or three pixels. 

3 3D and 4D TEXTURES 

We worked with both 3D and 4D procedural textures as input for advection of texture coordinates by 
the flow field. In this section, we describe how each sphere in the texture is represented, how the spheres 
are spaced, and how they change over time. We also describe the texture advection process. 

3.1 Spherical 3D and 4D Textures: 

Each 3D texture is populated with spheres. We initially use spheres to see how well the method 
works. We eventually, plan to replace the spheres with more cloud-like shapes. Currently, each sphere 
is assigned a texture value of 255 at the enter. This value linearly drops off to 127 as we go out along the 
radius. Outside the sphere, the texture values are set to 0. This allows for fuzzy spheres. To introduce 
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some variation, one can conceivably add noise to the values within the sphere, or change the range of 
values defining the sphere, etc. 

The spheres populating the 3D texture are randomly distributed throughout the volume using a poisson 
distribution. That is, we specify a poisson radius in addition to the sphere radius. The first frames of 
Figures 1 and 2 show the initial distribution of spheres in a 3D and 4D spherical textures respectively. 

In our current implementation, we statically generate the 3D texture volumes. To allow the introduc- 
tion of new spherical textures across inflow boundaries, we set the texture dimension to be twice the flow 
data volume along each dimension. This works for relatively short streamlines or animations. However, 
we will eventually run out of sphere textures to feed into the flow volume using this strategy. 

In our earlier work with 2D flow visualizations using this texture advection technique, we noticed 
that in some flows, the advected textures may get crowded in some regions of the flow, while in other 
regions, there is a sparsity of textures. One solution is to add an extra dimension to the texture space 
(3D textures for 2D flow fields) so that new textures may pop out along the new dimension to avoid the 
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sparsity problem. Likewise, old textures may disappear over time as they age in order to reduce the over 
crowding problem. So, in this work, we also consider 4D spherical textures. One can simply treat the 
4th dimension as time, where old spheres disappear over time, and new spheres appear over time. 

While this idea seem logical, it is not very practical due to the large storage requirement for 4D 
textures. Instead, we use additional frames of 3D textures as a bank of additional texture to draw from. 
In the examples shown in this paper, we use a bank of 4 3D texture frames, and cycle through these 4 
volumes as needed. The next section discuss how these textures are used. 

3.2 Advection of 3D and 4D Textures 

We use the VisTech library [1, 8] to perform backward integration of streamlines. It handled such 
tasks as point location so that we perform the advection in physical space and it finds the corresponding 
computation cell that contains the data. Since we are dealing with steady flow data sets, the integration 
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length of the streamlines also corresponds to the animation time. 

Here is the algorithm we used to advect a single 3D texture: 

1 . First, map the flow field to 3D texture space. 

2. For each 3D voxel, according to its physical position in the flow field, seed and initiate backward 
integration of a streamline of length ti. 

3. For each point along the streamline, obtain the corresponding texture coordinate and value (via 
trilinear interpolation), and assign it to the seed voxel. 

4. The streamline integration may terminate in one of three ways: 

(a) If the streamline has a length of zero (e.g. for voxels at the boundary or at source critical 
points), then a texture vaiue of zero is assigned to each frame in the animation for this voxel. 

(b) If the streamline reached a boundary or a source critical point after fnoving num steps, but 
still less than the total of n steps, then introduce n — num voxels from the texture volume. 

(c) If the streamline was entirely within the texture volume, then the corresponding texture val- 
ues along the streamline are assigned to the seed voxel in consecutive frames of the flow 
animation. 

5. Output advected texture volume as a time frame in the flow animation. 

Figure 1 shows a few frames from using 3D spherical textures for a synthetic flow field where the 
flow goes around a ring to produce a donut flow pattern. Note that no new spherical textures can be seen 
entering the flow volume. For that, we use 4D textures. The algorithm is essentially the same as above, 
except for step 4(b). 

If the streamline reached a boundary or a source critical point after moving num steps, but still less 
than the total of n steps, then introduce n — num voxels from the texture volume. If there is less than 
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n — num texture voxels available in the current 3D texture frame, then we advance to the next 3D texture 
frame. If we run out of texture frame, we simply recycle the 4 frames over again. Figure 2 illustrates the 
effect of 4D textures on the same vortical ring data set. 

4 VOLUME RENDERING 

We use direct volume rendering (DVR) to render the sequence of 3D advected textures that have 
been advected by the 3D flow field. From the sequence of rendered volumes, we collate them into 
an animation. Direct volume technique (DVR) is one of the most popular methods for visualizing 3D 
data sets. Some of the fundamental work can be found in Drebin et al. [3] and Levoy [7], More 
recent work, particularly on speeding up the algorithm, includes hardware texture mapping technique 
[2, 1 1] and shearwarp [6], For this project, we use a publicly available DVR software rendering package 
called UltraVis (www.hpl.hp.com/ultravis) from Hewlett Packard Laboratory. It is optimized for Intel 
Pentium III based PCs and achieves impressive interactive rendering rates on volume data of reasonable 
sizes. We chose this program because it is fast, freely available, does not depend on hardware and uses an 
accurate raycasting method. For animation purposes, we defined transfer functions so that the spherical 
texture patterns have fuzzy boundaries in order to simulate motion blur effects. We did not use any 
surface shading effects because they do not enhance the perception of the flow' considerably but rather 
exaggerate the artifacts on the images due to the low resolution of the data. In order to give some depth 
of field, we also defined semi-transparent materials for the empty space between the spherical texture 
patterns. 

5 RESULTS 

We tested our 3D and 4D texture map with 3D flow field data. In this section, we describe some of 
the animations we generated. 
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Figure 1 and 2 shows snapshot images from the animations of the vortex ring data. The data dimension 
is 64 3 and has vortex ring where the flow goes around from the top of the volume, throught the central 
axis and back out and over the top like a donut. In Figure 1, we use a single 3D texture volume. From 
the animation sequence, we observe that the spheres get stretched as they are advected over time. The 
stretching simulates motion blur effects, and can also be enhanced using different transfer functions 
during rendering. However, over-stretching would also hinder the desired effects. So, an appropriate 
balance must be struck. 

In Figure 2, we use multiple 3D texture volumes, we notice that some of stretched spheres disappear, 
while new spheres are introduced. This type of texture seem to enhance the dynamic nature of the flow, 

and also reduce the problem of over-crowding of textures and sparsity of textures in different regions of 
the flow over time. 

In Figure 3, we apply the same techniques to a CFD data set that describes flow over and around an 
airfoil. The dimension of this data set is 1 15 x 157 x 83. We also use multiple ^D texture volumes. For 
the wing geometry, we assigned a pre-arranged texture value so that the wing can be rendered together 

with the flow. This helps provide a frame of reference, but we still need to improve on smoothing the 
geometry. 

6 SUMMARY AND FUTURE WORK 

We presented some initial results on how one might extend texture advection to 3D flow volumes. 
We have found that straight forward extensions of 2D techniques is possible but may not be the most 
efficient in terms of storage requirement. The results are still quite early but do show some promise. 
There are a number of avenues that we are currently exploring to improve this work. This includes: (1) 
make sure that proposed technique is robust enough to handle other types of 3D critical points and/or 
curves, (2) explore other types of 3D and 4D textures. In particular, our goal is to produce cloud like 
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textures because most viewers can relate to their experience of watching clouds blown by the wind and 
inferring the wind pattern from their appearance. Rather than seeding or growing the cloud textures, we 
are also investigating an alternative subtractive process where we remove “matter” from a solid texture to 
produce cloud-like textures; (3) experiment with the rendering parameters to help produce a more fuzzy 
and wispy looking rendering of the advected cloud textures; (4) improve the space and time requirements 
of the algorithm; and (5) extend technique to time-varying 3D flow fields. 

ACKNOWLEDGMENTS 

We would like to thank Chris Henze for the synthetic 3D vortex ring flow field, Jennifer Dacles- 
Mariani and Greg Zilliac for the CFD wing data set, Suzana Djurcilov on alternative cloud-like 3D 
textures, and Gunter Knittel and HPL for making Ultra Vis available to the community. We would also 
like to thank the members of the Advanced Visualization and Interactive Systems laboratory at UC, Santa 
Cruz for their feedback and suggestions. This project is supported in part by N^SA grant NCC2-5281, 
LLNL Agreement No. B347879 under DOE Contract No. W-7405-ENG-48, NSF NPACI ACI-96 19020, 
and NSF ACI-990888 1 . 

References 

[1] Stephen T. Bryson, David Kenwright, and Michael Gerald- Yamasaki. FEL: The Field Encapsulation Li- 
brary. In R.D. Bergeron and A.E. Kaufman, editors, Proceedings of Visualization 96, pages 241-247. ACM, 
October 1996. 

[2] Brian Cabral, Nancy Cam, and Jim Foran. Accelerated volume rendering and tomographic reconstruction 
using texture mapping hardware. In Proceedings 1994 Symposium on Volume Visualization, pages 91—98, 
Washington, D.C., Oct 1994. IEEE/ACM. 

[3] R. A. Drebin, L. Carpenter, and P. Hanrahan. Volume rendering. In Proceedings of SIGGRAPH 88, pages 
65-74, August 1988. 


9 



[4] A.J.S. Hin and F.H. Post. Visualization of turbulent flow with particles. In Proceedings of Visualization 93, 
pages 46-52. IEEE, 1993. 

[5] David Kao and Alex Pang. On animating 2D velocity fields. In SPIE Visual Data Exploration and Analysis 
Vlll y 200 1 . to appear. 

[6] Philippe Lacroute and Marc Levoy. Fast volume rendering using a shear-warp factorization of the viewing 
transformation. In Proceedings of SIGGRAPH 94, pages 451-458, Orlando, FL, July 1994. 

[7] Marc Levoy. Display of surfaces from volume data. IEEE Computer Graphics and Applications, 8(5):29-37, 
May 1988. 

[8] P. Moran, C. Henze, and D. Ellsworth. FEL 2.2 user guide. Technical Report NAS-00-002, NASA, 2000. 

[9] Han-Wei Shen and David L. Kao. UFLIC: a line integral convolution algorithm for visualizing unsteady 
flows. In Proceedings of Visualization ’97, pages 317-322, 556. IEEE, October 1997. 

[10] Johan Stolk and Jarke J. van Wijk. Surface-particles for 3D flow visualization. In Proceedings Second 
Eurographics Workshop on Visualization in Scientific Computing, 1991. 

[11] Allen Van Gelder and Kwansik Kim. Direct volume rendering with shading via 3D textures. In ACM/IEEE 
Symposium on Volume Visualization, pages 22-30, San Francisco, CA, October 1996. 

[12] R. Wegenkittl, E. Groller, and W. Purgathofer. Animating flow fields: rendering of oriented line integral 
convolution. In Computer Animation '97, pages 1 5-2 1 . IEEE Computer Society Press, June 1997. 


10 




frame 17 frame 30 


Figure 1. Images from the animation of vortex data with single 3 dimensional texture map. 
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frame 27 frame 34 


Figure 2. Images from the animation of vortex data with 4 dimensional texture map. 
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Figure 3. Images from the animation of a 
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frame 30 

FD wing data using four 3D texture volumes. 
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